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Abstract 



An analytical model is presented for a rectangular lattice of isotropic scatterers 
with electric and magnetic resonances. Each isotropic scatterer is formed by 
putting appropriately 6 f2-shaped perfectly conducting particles on the faces 
of a cubic unit cell. A self-consistent dispersion equation is derived and then 
used to calculate correctly the effective permittivity and permeability in the 
frequency band where the lattice can be homogenized. The frequency range in 
which both the effective permittivity and permeability are negative corresponds 
to the mini-band of backward waves within the resonant band of the individual 
isotropic scatterer. 

PACS: 41.20.Jb, 78.20.Ci, 42.25.Bs, 42.70.Qs. 
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1 Introduction 



In 1967 Veselago considered an isotropic continuous medium with both negative 
permittivity and negative permeabihty and proved that a negative refraction 
would occur when a monochromatic wave impinges on the surface of such a 
medium [1]. As a result he obtained theoretically a quasi- lens (the lens does not 
focus parallel rays) from a slab of such a hypothetic material. In the Veselago 
theory this quasi-lens must form an image of any source without distortion (since 
each interface is a parallel plate). The phase of the radiated field at the focal 
point should be the same as the phase of the source dipole. These results are the 
consequences of backward waves in such a medium. The phase velocity and the 
Pointing vector are in opposite direction for a backward wave. In 2000 Pendry 
claimed that the Veselago quasi-lens must possess an extraordinary property 
to amplify and focus the evanescent waves in the spatial spectrum of a point 
source and therefore this lens must be perfect [2]. Since then much attention 
has been attracted to this research area both theoretically and experimentally 
(see e.g. [3] — [7]). Such negative media do not exist in nature but have already 
got different names in the literature such as backward- wave media, the Veselago 
media, left-handed media and negative mcta-materials. Since it is not clear yet 
how to make such media in the optical range of frequencies, the experimental ef- 
forts have been concentrated on creating negative meta-materials at microwave 
frequencies. A structure allowing the propagation of backward waves within a 
certain frequency band has been suggested in [4] and then experimentally stud- 
ied in [3], where the negative refraction phenomenon was primarily observed. 
However, the structure of [4] (a combination of two lattices: a lattice of in- 
finitely long parallel wires and a lattice of the so-called split-ring resonators) is 
not isotropic. This structure can be described with negative permittivity and 
permeability only if the propagation direction is orthogonal to the axes of wires 
[5]. 

To the best of our knowledge, the only isotropic medium with possible neg- 
ative permittivity and permeability (within a certain frequency band) that has 
been reported is formed by many isotropic cubic cells of Q particles [6]. Each 
isotropic cubic cell is made by putting 6 Jl-shaped perfectly conducting particles 
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on its faces (as shown in Fig. 1). Fig. 1 shows how to make an isotropic unit 
cell by putting successively 3 pairs of f2-particles on the opposite faces of a cubic 
unit cell in an appropriate way. Then each cubic cell can be described approx- 
imately as an isotropic resonant scatterer. The frequencies at which both the 
effective permittivity and permeability of such a composite medium are negative 
arc within the resonant band of each individual isotropic scatterer. Note that 
conducting fi-particles were first suggested in [8] for creating the bianisotropic 
composites. 

In the present paper we consider a negative meta-material formed by a rect- 
angular lattice of isotropic cubic unit cells of Q particles (as shown in Fig. 1). 
We introduce an analytical model to describe the dispersion properties of the 
lattice and calculate correctly the effective permittivity eeff{u)) and permeabil- 
ity i^effi^^) in the frequency band where the lattice can be homogenized. We 
compare the results obtained by the present model with the previous results 
obtained by the Maxwell Garnett model, and show that the latter can not give 
reliable values for either the real parts or the imaginary parts of the effective 
material parameters /Xe// and Ce//. 

Prom the negative dispersion of the lattice at low frequencies we obtain the 
frequency range of backward waves within the resonant band of the isotropic 
scatterers. In this frequency range the phase velocity is in an opposite direction 
of the group velocity. Our calculations show that this backward wave range is 
exactly the same as the frequency range within which both Cgff and /Xg// are 
negative. 

The case when the lattice is formed by many parallel bianisotropic f7-particlcs 
is also considered in an appendix. It is shown that the bianisotropy of the lattice 
particles reduces the possibility for the existence of backward waves. 

2 Polarizability of an isotropic cubic unit cell 

For bianisotropic particles (a 0-particle is a special case), one has the following 
four dyadic polarizability matrices relating the electric dipole moment p and 
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the magnetic dipole moment m with the local electric and magnetic fields: 

p = iee•El°=+lem•H'°^ (1) 

m = i„e-Ei°^+i„„-H'°^ (2) 

For reciprocal particles, one has a^e = —o-me- 

For a cubic unit cell of O particles as shown in Fig. 1, the situation is 
quite different. The bianisotropy cancels out since e.g. each pair of opposite 
O particles shown in Fig. 1 have equal absolute value but opposite sign for 
aem- The polarizabilities ciee and 'amm of a cubic unit cell of 6 f2-particles 
(see Fig. 1) are the sums of the polarizabilities a^; and for individual 

particles [6] . The influence of the mutual coupling of the particles is small and 
the resonant frequency for such a cubic unit cell of 6 fi-particlcs has only a 
small shift from that for an isolated il-particlc. The resonant excitation by an 
electric field directed along the arms of il particles is mainly due to the presence 
of the arms. The resonant behavior retains for the total electric polarization of 
the two opposite omega particles. Within the resonant band, the quasi-static 
polarizability component ag|^ (which describes the response of the particle to 
an electric field normal to the direction of the arms; see Fig. 1(a)) of the O 
particle is small as compared to a|^^. Such a cubic unit cell can be therefore 
considered as an isotropic scatterer and thus one has = a^^ = aH = Uge 
and = = = amm [6] . One can easily show that [6] 

flee ~ 2(Xgg , flmm — "^Oimm- i"^) 

The expressions for the parameters afg^ and al^^ have been given in [9]. In 
the present paper we consider a material formed by a rectangular lattice of such 
scatterers (see Fig. 1(b)). 

The electric and magnetic polarizabilities for a cubic cell of f2-particles are 
calculated (cf. [6]) and shown in Fig. 2 as a numerical example. The geometric 
parameters for the 0-particles are chosen as r = 1.5 mm, w = 0.4 mm, h = 0.2 
mm and I = 2 mm (see Fig. 1(a)). The size of the cubic unit cell is 4 mm. 
The relative permittivity for the background medium is chosen as = 1.5 (the 
permeability of the background medium is assumed be the same as the one for 
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vacuum in the present paper). Resonances of Uge and amm (frequencies at which 
the real parts of these parameters become zero) occur at the frequency of 8.05 
GHz and 8.14 GHz, respectively. We have chosen the time dependence e^"*, and 
thus the imaginary parts of these parameters are negative. Below we consider 
a material formed by a rectangular lattice of such isotropic scatterers (see Fig. 
1(b)). 

3 Wrong results predicted by the Maxwell Gar- 
nett model for the regular lattice 

In [6], the well-known Maxwell Garnett model was used and the following ex- 
pressions for the material parameters were obtained, 

1 f NUee N'^aeeamm\ , 
^eff = £6 + ^ > (4) 

where 

3eo 9eoe5/xo 
In the Maxwell Garnett model, one only needs to know the scatterer polar- 
izabilities and the density N of the scatterers. Thus the obtained result for a 
regular lattice of scatterers will be the same as that for a random distribution 
of scatterers if the density of the scatterers is the same. For the case considered 
in the present paper (i.e., a rectangular lattice of isotropic cubic unit cells as 
shown in Fig. 1(b) with periods dx,dy,dz along the Cartesian axes), one has 
N = l/dxdydz- The frequency dependencies of eeff and /ie// predicted by the 
Maxwell Garnett model are shown in Fig. 3 when d^ = dy = dz — d — 8 mm 
(i.e. N = (1/8^) • 10~^ m~^). From this figure one sees that both Re(/;ieff) and 
Re(eeff) are negative within a certain frequency band (8.34—8.62 GHz). How- 
ever, within this frequency band the imaginary parts of fXeS and Cgs have high 
values (about —4 ~ —6). The radiation losses are so large that the obtained 
negative effective permittivity Cg/ / and permeability /ig/ / have no practical im- 
portance for applications. 



6 



The above results predicted by the Maxwell Garnett model (in [6] for a 
random arrangement of cubic unit cells) are not correct for the present case of 
the regular lattice. Within the resonant band (7.6-8.6 GHz) of each il-particle, 
the wavelength (A = 28 — 33 mm) in the background medium (et = 1.5) is 
much larger than the distance d (8 mm) between the isotropic ciibic scatterers, 
and this distance is larger than the size (4 mm) of each scattcrcr. From this 
point of view one may expect that the homogcnization is possible. However, 
the results predicted by the Maxwell Garnett model arc inconsistent near the 
resonance of the scatterer. The Maxwell Garnett model is based on the Clausius- 
Mossotti relations for the electric and magnetic fields. These relations for cubic 
lattices of static dipoles have been found as good approximations in the works 
of Sivukhin [10] (for an infinite lattice and for a half-space) and McPhedran [11] 
(for arrays of finite sizes). The Clausius-Mossotti relations were also confirmed 
to be accurate enough for cubic lattices of dipoles in a time-harmonic case [12] 
under the condition that the wavelength is large compared to the lattice period 
d. However, the resonant case is not considered in these works. In the resonant 
case, the Maxwell Garnett model may give very high absolute values for the 
real parts of Cg// and i^eff within the resonant band, and this would lead to 
a dramatic shortening of the wavelength in the homogenized medium. The 
wavelength in the homogenized medium becomes comparable with the particle 
distance d and thus the Maxwell Garnett model becomes contradictory at some 
frequencies. 

In the present case (sec Fig. 3) the absolute values of real parts of e^ff 
and fJ^eff are so high within the narrow band 8.37-8.40 GHz where the Maxwell 
Garnett model becomes contradictory. Outside this narrow band |Rc(eejy )| and 
|Re(y^e//)| are not very high. However, the inconsistency of the Maxwell Garnett 
model for the present case of regular lattice is evident over the whole resonant 
band of f2-particles (8.1-8.5 GHz). In Fig. 3 one can see that the imaginary 
parts of Ceff and /Ug// are rather significant at those frequencies. On the other 
hand, we know that lm{ij.es) and Im(eeff) must be identically zeros physically if 
the lossless scatterers are arranged in a regular lattice. The cancelation of the 
radiation resistances of the scatterers (due to their electromagnetic interaction) 
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is obvious for an infinite regular array and thus there is no radiation loss at all 
(see e.g. [13]). 

In the present paper, we introduce another homogenization model which 
gives correctly zero value for the imaginary parts of Ce// and fieff- A com- 
parison in Section 5 will show that the real parts of Ce// and fJ^eff predicted 
by the Maxwell Garnett model are also wrong within the resonant band of the 
scatterers. 

4 Dispersion equation for the lattice of the isotropic 
scatterers 

Consider the rectangular lattice of the isotropic scatterers (shown in Fig. 1) with 
periods , dy along the x and y axes and dz — d along the z axis. In the present 
case, the electric dipole moment is orthogonal to the magnetic dipolc moment 
for each scatterer. Let p and m denote the dipole moments of a reference 
scatterer located at the origin. Thus one has 

p = aeeE^-, (6) 
m = a^^H'-. (7) 

We need to find the cigenniodcs for such a lattice of isotropic scatterers. Assume 
the wave propagates along the 2;-axis. Let (3 denote the propagation constant 
of an eigenmode. The first Brillouin zone is the interval < (3d < ty. Then 
p{nx,ny,nz) and m(n2;, n^, n^) (the electric and magnetic dipole moments of 
a scatterer located at a lattice node with coordinates x = nxdx,y = nydy,z = 
Uzd) can be expressed through p and m (for the reference scatterer) as 

P(nx,%,n,) = pe-^"^''^ (8) 
m(n^,%,n^) = me"^"-'''^. (9) 

Since eigenwavcs are linearly polarized in the present case, we can assume 
without loss of generality that p = pxo and m ~ myo (sec also Fig. 1). Then 
the local electric and magnetic fields arc directed along the x and y directions, 
respectively (the local field has the same polarization as the eigenmode). This 
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allows the following scalar expressions for the local field amplitudes in terms of 
p and m: 

E^"" = Ap+Bm, (10) 
H^°'= = Cm + Dp, (11) 

where A, B, C, D are the so-called interaction factors of the lattice, which depend 
only on the lattice geometry, the frequency co and the propagation factor /3 (they 

do not depend on the polarizabilities of the scatterers). From the reciprocity 
it follows that B = D. From the duality it follows that C = A/rj^, where 



V = V A*o/£b£o is the wave impedance of the background medium. 

From Eqs. (6), (7), (10) and (11), one obtains the following two important 
relations, 

p{l - tteeA) = aeeDm, (12) 

m{l - ^haf^) = ammDp. (13) 
Denote the dimensionless ratio r]p/m as a, then we have 



/ _!1 A 

p _ V^™™ V 



a = rj- = ^ . (14) 

This formula will be used in the next section to find the effective material 
parameters of the lattice at low frequencies. 
From Eqs. (12) and (13) one obtains 

(—-A)(^-A)=r,'D. (15) 

\ ^ee / \ ^mm / 

The relation 

lm(l-)= ^^^h!^ (16) 
\ flee / DTreo 

was first obtained in [13] for a particle re-radiating the light and later was 
reproduced as the consequence of the energy conservation for a lossless dipole 
scatterer (see e.g. [14]). For lossless magnetic scatterers, one has the following 
corresponding form, 



Im ( ^— ^ = ^^^^i^oetlio) ' 
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For the interaction factor A, an explicit approximation has been found in [14]: 

nu ( sin kd \ I? , ^ 

^ = ^n-^W^ — tta irA^i-a > 

zdxdy \ cos kd—cospd J oneoCb 
where k = w^eo/^oCb is the wavenumber in the background medium, and go is 
the real part of the dimensionless interaction factor of a 2D grid of dipoles with 
periods dx,dy. qo has the following explicit expression when dx = dy = a [14], 

1 / coskas . , A /,„s 
9o=oMz:— ' (18) 



2 \ kas 

where s w 1/1.4380 ~ 0.6954. Relation (17) is almost exact when d ^ dx,dy, 
and has only a small error when d = dx = dy [14] . 
Since 

1 \ [if 
\mA = Im ( = Im 



, a.ee 

the left-hand side of (15) is real. Thus wo have 



Re ( — - A Ve -A] = V'D. (19) 



Below we show that D is also real and we derive its explicit expression. From 
the definition of D, one has 

Hl"-{x = 0,y = 0,z = 0)=Y,Y.Yl Hyiuxdx, riydy, n,d) = Dp, (20) 

«i riy 

where Hy{nxdx, nydy, rizd) is the y-componcnt of the magnetic field at the origin 
produced by the x-polarized electric dipole with dipole moment 'p{nx, ny,nz) = 
p{nx, riy, nz)x.(). From the reciprocity we know that Hy{nxdx, Uydy, rizd) is equal 
to the field at the point {x = Uxdx, y = riydy, z = Uzd) produced by the reference 
dipole p. In Eq. (20) the summation is over all integers {nx,ny,nz) from 
— oo to oo except Ux — Uy — riz = 0. Here we use the following plane-wave 
representation for the magnetic field produced by an electric dipole p = pxo 
located at the origin [16]: 

Hy{x,y,z) = ^-^{z)p J J e-^^'''^+''yy^-''^'W'''-oi-oi^dqxdqy, (21) 

— oo — oo 

where *(2) = +1 for ^ > 0, ^{z) = -1 for z < and ^(z) = for z = 0. 
Substituting Eqs. (21) and (8) into definition (20) and changing the order of 
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the summation over nz and the integration over and qy, one obtains 

oo oo oc 
n«_oo-oo "»=1 

Summarizing the geometrical series and using known exphcit formulas for the 
integration over q^ and qy, one can obtain the following final result for D, 



ria, riy ^.Qg 



\lk'-{'ity-{^yd-cospd 



Numerical calculations have shown that the contribution from the terms 
with Tlx, riy = ±1, ±2, . . . becomes significant only if d^^dy d. These terms 
correspond to the effect of the evanescent Floquct spatial harmonics generated 
by the 2D grids of dipolcs with = ±1, ±2 . . . (the reference dipolc is located 
in the 2D grid with = 0) . If we consider the case d^ = dy = a < d, then we 
can neglect these terms and take into account only the plane-wave interaction 
between the 2D grids of dipoles. For such a case, one has the following real- 
valued expression for D, 

2dxdy cos kd — cos pd 
Note that relation (17) was derived in [14] under the same condition. Formu- 
las (17) and (23) are both approximate but self-consistent. Therefore, these 
relations lead to the real- valued dispersion equation (19). 
Substituting (17) and (23) into (19), one obtains 

cos^ Pd{l + 7172) — cos l3d (27172 cos kd — (71 + 72) sin kd) 

— (1 — 7172) cos^ kd — (71 + 72) sin fcd cos fed = 0, (24) 

where 



Re ^dxdy _ Re 2dxdy 



71 = , 72 



The dispersion relation (24) has two roots for a fixed frequency. One of the 
roots satisfies the following condition (obtained from Eqs. (12)and (13)) 



V amnv V ) _ 

^ " (--A' 

\aee ) 
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This root is /3 that we show in our dispersion plots below. The other root does 
not satisfy the above condition and thus is spurious. 

Fig. 4 shows the dispersion for the same cubic lattice of isotropic scatterers 
as used for Fig. 3. From this figure one can see that there is a band (between 8.14 
and 8.37 GHz) in which the group velocity is negative. This is the backward- 
wave band. Since the Foynting vector in a lossless medium must be directed 
along the group velocity, this band corresponds to the case when the energy 
transports oppositely to the phase. Of course, this situation would be trivial 
for some high-frequency dispersion branches. For a lattice of dipoles, negative 
dispersion is inherent for every even dispersion branch (7r(2n-|-l) < kd < TT2{n+ 

1) 7r, where n = 0,1,2...). However, only in the first frequency zone (where 
< kd < n) the lattice period is smaller than A/ 2 (half of the wavelength in 
the background medium) and thus the homogenization of the lattice is possible. 
Only if the lattice of scatterers can be homogenized, the material parameters 
Cef f and f can be introduced and the lattice of scatterers can be interpreted 
as a continuous medium. It the Veselago theory [1] the backward wave in 
a continuous medium should correspond to a case when both Ce// and /Xe// 
are negative. The negative dispersion in the homogenized structure implies 
the backward wave and therefore our eg// and /ig// should be negative if the 
Veselago theory is applicable for dispersive media. 

Fig. 5 gives an enlarged view of Fig. 4 for the frequency dependence of 
Re/3 and lm/3 in the frequency range 7-10 GHz. From this figure one sees that 
there are two stopbands within the resonant band of each scattcrcr (cf. Fig. 

2) . The lower one (stopband 1) is very narrow (8.06-8.14 GHz), and the higher 
(stopband 2) is wider (8.37-8.52 GHz). Between these two stopbands there is the 
band of backward waves. Stopband 2 is a conventional lattice stopband which 
has Re/3 = 0. Stopband 1 corresponds to the so-called complex mode which is 
known in the theory for lattices of infinite wires (see e.g. [17], [18]). This complex 
mode is decaying and the real part of the propagation factor is identically equal 
to n/d. In our case, the relation Re/3 = ir/d reflects the fact that the directions 
of the dipole moments of the scatterers are alternating along the propagation 
axis (i.e., two adjacent isotropic scatterers have opposite polarizations). 
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Substituting relations (17) and (23) into Eq. (14), one obtains the following 
explicit formula for the parameter a, 

go — -%T^Re ( ) (cos fed — cos /3(i) + sin fcrf 

a{uj) = ^ ^ ""^j ^. . (25) 

sm jjd 

Here amm is a function of w for given geometrical parameters of scatterers and 
given Cfc. Since the frequency dependence of is known from the solution 
of the dispersion equation (24), the frequency dependence of a{(jj) can be cal- 
culated from the above equation. In the next section we will use this frequency 
dependence of a{w) to find the explicit expressions for the eflfective material 
parameters ee//(w) and ^effico) of the homogenized lattice. 

Our calculations show that a is also a crucial parameter for the frequency 
bounds of the backward-wave range. The central frequencies of both stopbands 
are frequencies at which a{u) changes the sign. It is positive below 8.1 GHz, 
negative in 8.1-8.45 GHz (backward-wave range) and positive again above 8.45 
GHz. Both stopbands are centered by the 2 zero-points of a{u)). Of course, 
both stopbands are within the resonant band of each isotropic scatterer. 

5 Homogenization of the lattice of the isotropic 
scatterers 

Consider the frequency band < kd < it within which homogenization is pos- 
sible and effective parameters eg// and Heff can be, perhaps, introduced. To 
determine the frequency dependence of the two parameters eeff{Lv) and fj.eff{i^), 
we needs to find two relations (involving the propagation constant f3{(jj)) between 
these two parameters. The first relation can be obtained by fitting the value 
P/k with the eflfective refraction index rie// = y/CefffJ-eff, i-e., one has 

= ^/ 

Note that k is proportional to co (i.e., k = ui^/eojloeb) and /? = is given by 
the dispersion curve calculated in the previous section. 

To find another relation, we define the following averaged (over a cubic unit 



13 



cell) electric and magnetic polarizations per unit volume: 

where V = (in this section we assume that the lattice periods are the same 
in all the three directions, i.e., dx — dy — = d; otherwise the lattice will be 
anisotropic). In the same way we can introduce the following averaged electric 
and magnetic fields, 



<E>=^jEdV, <H>=^j- 



HdV. 

V V 

Since the material considered here is isotropic, its wave impedance Z can be 
expressed as the ratio < E > / < H >, i.e., 

Z. (27) 



<H> 



If we use the conventional expression Z ~ y ; the square root can be 

chosen either positive (in the range of forward waves) or negative (in the range 
of backward waves; the ratio fieff/^eff is always positive). For a backward wave 
in a continuous isotropic medium, the vectors < E >= Exq and < H >== -ffyo 
form a left-hand triad with the wave vector K — (3zo [1], and thus the wave 
impedance should be negative. Substituting (26) into Z = ^/^j^j^, one obtains 



(28) 



17^ P 



The above expression is more convenient for use since it gives automatically a 
negative impedance for a backward wave (when eg// < 0). 

From the definitions for the effective material parameters, one has 

eoee// <E>= eocb < E > +P, MoMe// < H >= < H > +M. 

Using Eqs. (27) and (28) and taking into account the above relations, one can 
express the dimensionless parameter a = r/P/M (determined by Eq. (25)) as 

eojceff-eb) <E> ^ 
" "^M '^Mo(Me//-l)<^?> 

^ eo(ee// - eb) ^ ^ {ceff - ^b) P ,^g. 

Mo(Me// - 1) V^ifJ-eff - 1) kCeff ' 
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The above equation gives the second required relation between the two effective 
parameters. 

Finally, from the two relations (26) and (29) we obtain the following expres- 
sion for determining ee//(w) from the dispersion relation (3{w), 

<^eff['^) - ^ p -■ (>}U) 

Note that k = u^eono^b and a = a{u)) is given by formula (25). After ee//(w) 
is calculated, iJ,eff{u)) is then calculated from Eq. (26). 

Fig. 6 shows the effective permittivity and permeability for the same lattice 
of isotropic scatterers as used before in Figs. 3—5. From this figure one sees 
that both Ceff and fieff are negative in the band of backward waves. This is 
consistent to the Veselago theory. 

In any of the two stopbands, the value /3 contains non-zero imaginary part 
and thus Eqs. (30) and (26) give complex values for Cg/f and fXeff at these 
frequencies. Since the structure is lossless, eg// and Heff must be real at the 
frequencies for which the homogenization is possible. Therefore, we can conclude 
that the homogenization can not be performed within the stopbands. That's 
why in Fig. 6 no value is given for eeff or fJ,eff in the stopbands. From Fig. 
6 one sees that outside the two stopbands the absolute values of e^ff and /ie// 
are not very high. The wavelength in the homogenized medium (which is equal 
to 2tt / (3) is not dramatically shortened within the passbands (in which [3d < 
tt). This indicates that the present homogenization model is self-consistent 
(unlike expressions (4) and (5) derived from the Maxwell Garnett model). By 
comparison Fig. 6 with Fig. 3, we can conclude that for a lattice of resonant 
particles the Maxwell Garnett model can not give reliable values for either the 
real parts or the imaginary parts for the effective material parameters (/Ue//) and 
(ce//). The Maxwell Garnett model can only predicts qualitatively the location 
of the frequency band where the effective permittivity and permeability have 
negative real parts. 
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6 Conclusion 



In the present paper, we have presented an analytical model for an isotropic 
negative meta-material, which is formed by a rectangular lattice of isotropic 
scatterers with electric and magnetic resonances. Each isotropic scatterer is 
formed by putting 3 pairs of f2-particles on the opposite faces of a cubic unit cell 
in an appropriate way. We have derived a self- consistent dispersion equation and 
studied the dispersion properties of the lattice. The obtained dispersion curves 
are then used to calculate correctly the effective permittivity and permeability 
in the frequency band where the lattice can be homogenized. It is interesting 
to see that the dispersion curves agree well with the Veselago theory which 
predicts backward waves when both permittivity and permeability are negative. 
The Veselago theory is non-complete (it was only for monochromatic waves and 
did not take into account the dispersion of the mcdiTim) and some questions are 
left open. It has been shown in the present paper that the frequency range in 
which both the effective permittivity and permeability are negative corresponds 
to the mini-band of negative dispersion of the lattice at low frequencies within 
the resonant band of the individual isotropic scatterer. It has been shown that 
the frequency range of backward waves is determined by the edges of the two 
stopbands at low frequencies. We have compared the results obtained by the 
present model with the previous results obtained by the Maxwell Garnett model, 
and shown that the latter can not give reliable values for either the real part 
or the imaginary part for any effective material parameter {i^eff or Ce//) of the 
lattice. For lattices of resonant scatterers, the Maxwell Garnett model can only 
predict approximately the location of the frequency band where the effective 
permittivity and permeability have negative real parts. 

The case when the lattice is formed by parallel bianisotropic fi-particles has 
also been studied (in the appendix). It has been shown that the bianisotropy 
of the lattice particles suppresses those terms in the dispersion equation that 
can allow a backward wave propagation and thus reduces the possibility for the 
existence of backward waves. Therefore, the present formulas for calculating 
correctly eg// and /^e// for the negative meta-material from resonant isotropic 
scatterers (the only isotropic negative meta-material that we know) is of impor- 
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tance. 
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Appendix: The case of a lattice of bianisotropic 
particles 

In this appendix, we present briefly our studies for the case when the lattice is 
formed by an inflnite number of parallel Q-particles. We also explain why the 
magnetoelectric coupling of these bianisotropic particles reduce the possibility 
for the existence of backward waves. 

We use a simplified wire-and-loop model for each Q-particle, in which we 
assume that the current flowing around the loop is uniform (i.e., the electric 
polarization of the loop is neglected). 

If the particle is excited by a a;-directed electric fleld (see Fig. 1), then 

the current in the loop and the straight wire arms will be equal to the induced 

voltage in the straight wire portion divided by the impedance for the series 

connection of the loop impedance Zi and the straight wire impedance Zyj, i.e., 

^_ £ _ E^^^l 
Zi + Z^ Zi + Zw 

Thus one obtains the following electric and magnetoelectric polarizabilities, 

_ II _ . P 
- jujEioc - ^Zi + Zj 

_ IhqS _ fXpSl 
ELoc - Zi + Z^' 

where S = 7rr^ is the loop area. 

Now let the particle be excited by a y-directed local magnetic fleld H''"'^. An 
electromotive force (i.e., voltage £) is induced in the loop, and one has 

Zl + Zy^ Zl + Zy^ 
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Thus, one has electromagnetic and magnetic polarizabilities 

_ II ^J'ols 



^ IiiqS _ jujiiqS 



Prom the above expressions for the polarizabihties, the following identities 
can be obtained: 

(A2) 

The second identity is due to the reciprocity property of the particle. For such 
a bianisotropic particle (local electric field is x-polarized and magnetic field is 
y— polarized), we have (cf. Eqs. (6) and (7)) 

p = aeeE^°'= + aemH^'"', (A3) 
m = a„e^^°' + amm^?'°". (A4) 

In the rectangular lattice of parallel fJ-particles, the electric dipole moment of 
each particle is directed along the x-axis and the magnetic dipole moment is 
directed along the y-axis when the eigenmode with E = i^xo and H — Hyo 
propagates along the axis (the local field has the same polarization as the 
eigenmode). Substituting relations (10) and (11) into Eqs. (A3) and (A4), one 
obtains 

p{l - UeeA - GemD) = m{aeeD + aem^), (A5) 
m{l - amm^ - CLmeD) = p{ammD + ttemA). (A6) 

The above two equations give 

A A 

(1 (leeA (lgrn-D){l (lmm~o (^meD) = {(IggD + (lem~o)(flmm-C' + ClemA). 

Since all terms containing D cancel out due to identities (Al) and (A2), one 
obtains 



1 ^eeA ^mm ^ — Oi 
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which can be rewritten in the following form 

1 



A. (A7) 



For a bianisotropic particle, Eq. (16) is generalized to the following relation 
[19], 

Im ( L\ = ^^f^. (A8) 

\aee + ^J 67reoef, ' 

Substituting Eq. (17) into Eq. (A7) and using Eq. (A8), one obtains the 
following dispersion relation, 

r, 1 , , sin /3d . , 

cos pa = cos ka H — — , (A9) 

where 

24(iyRe ( - qo 
X= ^ 



UJ1] 

Eq. (A9) has the same form as the classical dispersion equation for a periodically 
loaded line [20]), which is also known as the Kronig-Pcnnic equation in the solid- 
state physics. In the known cases described by this equation (see e.g. [20]) the 
negative dispersion exists only within the even zones of frequencies {2nn < kd < 
7r(2n + 1), n = 1,2,3, ...). We did not find any physically realizable O-particle 
for which the above dispersion relation would give a negative dispersion within 
the first frequency zone {0 < kd < it). 

In the above analysis, the model for the O-particle is rough (since the electric 
polarization of the loop is neglected) and identity (Al) is approximate. We have 
studied the lattice of parallel 0-particles using a more accurate model (as used 
in [9]). This accurate model leads to another dispersion equation (different 
from Eq. (A9)) which is quadratic with respect to cos (3d and contains the 
electromagnetic interaction factor D. Again, we did not find any physically 
realizable f2-particle for which this dispersion relation would give a negative 
dispersion within the first frequency zone {0 < kd < n). The error in Eq. (Al) 
for real Jl— particles is quite small and the terms containing D in this dispersion 
equation arc relatively small. Neglecting these small terms and removing the 
spurious root, we obtained the same result as what Eq. (A9) leads to. 
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Therefore, we can conclude that the bianisotropy of the lattice particles sup- 
presses those terms in the dispersion equation that can allow a backward wave 
propagation at low frequencies. The possibility for the existence of backward 
waves is thus reduced (it was reported with another reason in [21] that the 
bianisotropy is not good for the existence of backward waves). 
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Figure captions 




Figure 1: W (b) (a) Geom- 

etry for a cubic unit cell of rj-particles. For graphic clearness, only 1 pair 
of ri-shaped perfectly conducting particles are shown on 2 opposite faces 
of the cubic unit cell, (b) A rectangular lattice of isotropic cubic unit cells 
of f2-particles. 




Figure 2: Left: The 

real (solid lines) and imaginary (dashed lines) parts of the electric polar- 
izability (left figure) and the magnetic polarizability (right figure) for an 
isotropic scatterer (a cubic unit cell of il-particles) . 
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7 7.5 8 8.5 9 9.5 10 

Figure 3: frequency, ghz The fre- 

quency dependencies of ee// and /ie// predicted (incorrectly for the regu- 
lar lattice) by the Maxwell Garnett model for the negative meta-material 
formed by a lattice of isotropic cubic unit cells of 0-particles (as shown in 
Fig. 1(b)). Re(eeff) (solid line), Re(/Xeff) (dashed line), Im(eeff) (circles), 
and Im(/Zeff) (crosses) are shown. 




Figure 4: P*jt The dis- 

persion (indicated by the circles) for the same lattice of isotropic resonant 
scatterers as used for Fig. 3. The thin lines indicates the dispersion of the 
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homogenous background medium. 




Figure 5: frequency, ghz An enlarged 

view of Fig. 4 for the frequency dependence of Re/3 (crosses) and Im/3 
(circles) near the resonance. The thin line indicates the dispersion of the 
homogenous background medium. 




Stopbands 




Backward wave band 



7 7.6 8 8.5 9 9.5 10 

Figure 6: frequency, ghz The fre- 

quency dependencies of the effective permittivity eg// (circles) and the 
effective permeability fieff (crosses) for the same lattice of isotropic reso- 
nant scatterers as considered in Figs. 3—5. 
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